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We consider a two-level system such as a two-level atom, interacting with a cavity field mode in the 
rotating wave approximation, when the atomic transition frequency or the field mode frequency is 
periodically driven in time. We show that in both cases, for an appropriate choice of the modulation 
parameters, the state amplitudes in a generic n-excitation subspace obey the same equations of 
motion that can be obtained from a static non-Hermitian Jaynes-Cummings Hamiltonian with 
VT symmetry, that is with an imaginary coupling constant. This gives further support to recent 
results showing the possible physical interest of VT symmetric non-Hermitian Hamiltonians. We 
also generalize the well-known diagonalization of the Jaynes-Cummings Hamiltonian to the non- 
Hermitian case in terms of pseudo-bosons and pseudo-fermions, and discuss relevant mathematical 
and physical aspects. 
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I. INTRODUCTION 

Non-Hermitian Hamiltonians have been frequently 
used in physics to describe dissipative effects in concrete 
systems in many different fields of physics, in particu¬ 
lar in quantum optics 0 , 0 , or for describing decaying 
states 0 . Recently, there has been a renewed interest in 
non-Hermitian Hamiltonians after it has been shown that 
non-Hermitian Hamiltonians with VT (parity-time) sym¬ 
metry can have a real eigenvalues [J-|6|. The same hap¬ 
pens for non-Hermitian but pseudo-symmetric Hamilto¬ 
nians, Q, where the VT symmetry is replaced by a more 
abstract condition. In a recent paper, it has been shown 
that physical systems such as optical lattices with some 
external modulation of their parameters can be described 
by a Lee-Friedrichs Hamiltonian in the ghost regime, that 
is when the coupling constant is imaginary and there¬ 
fore the Hamiltonian is not Hermitian]sf, 9j. In this pa¬ 
per we consider another physical system, of interest in 
quantum optics, with a time dependence of a physical 
parameter that can be simulated by a time-independent 
non-Hermitian Hamiltonian with an imaginary coupling 
constant. 

We consider a two-level system, a two-level atom for 
example, interacting with a boson field, for example a sin¬ 
gle cavity mode of the electromagnetic field, when one pa¬ 
rameter of the system is periodically modulated in time. 
Specifically, we investigate two cases: a modulation of the 
transition frequency of the two-level system and a modu¬ 
lation of the frequency of the cavity mode. Cavities with 
a conducting wall subjected to a periodic motion and 
optically modulated cavities have been extensively con¬ 
sidered in the literature, in particular in the framework 
of the dynamical Casimir effect EMU- We consider 
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the interaction of the two-level atom with a single cavity 
mode in the rotating wave approximation, and thus the 
system is described by a Jaynes-Cummings(JC) Hamilto¬ 
nian (T2I [Tit , with one time-dependent parameter (atomic 
transition frequency or cavity field frequency). Our sys¬ 
tem is thus described by a Hermitian time-dependent JC 
Hamiltonian, and we obtain the differential equations for 
the state amplitudes in each subspace characterized by 
a well-defined number of total (atom plus field) excita¬ 
tions. This is possible because our Hamiltonian com¬ 
mutes with the total excitation-number operator. We 
show that, under an appropriate choice of the param¬ 
eters characterizing the system’s modulation and after 
a time-average, these equations are identical to those 
that can be obtained by a static JC Hamiltonian with 
an imaginary coupling constant, that is a non-Hermitian 
JC Hamiltonian. This Hamiltonian, however, still has 
a VT symmetry. This result gives a physical meaning 
and importance in considering a non-Hermitian Hamil¬ 
tonian for a matter-field interacting system, specifically 
a fermion field linearly coupled to a boson field. This 
is an extension of an analogous recent result considering 
a modulated Lee-Friedrichs Hamiltonian (that is in the 
one-excitation subspace), where light transport in a semi¬ 
infinite waveguide lattice with a nonlinear modulation 
was considered and equivalence to a non-Hermitian Lee- 
Friedrichs Hamiltonian in the ghost regime was shown 
0. Due to the physical relevance of a non-Hermitian 
JC Hamiltonian, we then generalize the well-known di¬ 
agonalization of the JC Hamiltonian by a unitary trans¬ 
formation to the non-Hermitian case in terms of pseudo¬ 
boson [Util and pseudo-fermion operators E2 , which 
have been recently introduced for fermionic and bosonic 
systems separately. These pseudo-Hermitian operators 
arise from a special deformation of canonical bosonic and 
fermonic commutation relations 0E3. The pseudo- 
Hermitian Hamiltonian we obtain is physically equiva¬ 
lent to the non Hermitian VT symmetric Hamiltonian 
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and can be diagonalized exactly, generalizing the proce¬ 
dure used in the Hermitian case. Thus, our results extend 
the definition of pseudo-bosons and pseudo-fermions to 
the case of a fermion-boson interacting system such as 
the Jaynes-Cummings model, allowing us a mathemati¬ 
cally rigorous treatment of our non-Hermitian interaction 
Hamiltonian. The relevance and application of the for¬ 
malism of the pseudo-Hermitian operators in the frame¬ 
work of open quantum system or in quantum optics has 
been recently discussed in the literature (see for example 
fl8{ and reference therein) 

This paper is organized as follows. In Sec. [TT]we intro¬ 
duce our VT symmetric Jaynes-Cummings model with 
a modulation of the atomic or field frequency. In Sec. 
1III1 we consider the equations of motion for the state of 
the interacting system in each subspace characterized by 
a well-defined number of excitations, and show that, for 
appropriate choices of the system parameters, they are 
equivalent to a static non-Hermitian JC Hamiltonian. In 
Sec. El motivated by our previous analysis, we diago¬ 
nalize the non-Hermitian JC Hamiltonian through an ap¬ 
propriate transformation operator, generalizing the well- 
known diagonalization for the Hermitian case, in terms 
of pseudo-boson and pseudo-fermion operators, and dis¬ 
cuss some mathematical and physical aspects. Sec. [V] is 
finally devoted to our concluding remarks. 


II. THE MODEL 

A two-level atom interacting with a single mode of the 
electromagnetic field in the rotating wave approximation 
(RWA), can be described by the Jaynes-Cummings model 
[12H14I ] . This Hamiltonian is frequently used in quantum 
optics, in particular for studying non-perturbative effects 
in strong coupling regimes [19], [20]. The relative Hamil¬ 
tonian can be written as 

H = H 0 + V = fiioo ^ + hwala + g s (atf + a)b) , 

(1) 

where the annihilation and creation operators b ', b re¬ 
fer to the two-level system with transition frequency wo 
and satisfy standard fermionic anti-commutation rela¬ 
tions, a t , a refer to the field mode with frequency ui and 
satisfy standard bosonic commutation relations, and g s 
is a coupling constant that in this moment we assume 
to be real, so that the Hamiltonian © is clearly Hermi¬ 
tian. The JC Hamiltonian is the simplest Hamiltonian 
model describing the coupling between a fermionic and a 
bosonic field, and it can be diagonalized exactly. In fact, 
it is straightforward to verify that it commutes with the 
operator giving the total number of excitations (atom 
plus field) 

A f = b^b + a) a . (2) 

Thus the Hamiltonian can be diagonalized independently 
in each subspace characterized by a well-defined number 


of excitations. In the n-excitation subspace, a basis is 
given by the two states | n, 0) and | n — 1,1), where 
the first element in the states indicates the number of 
excitations in the bosonic sector, and the second one that 
in the fermionic sector. The first state has n quanta in the 
field mode and the two-level system is in its ground state, 
while in the second state there are n— 1 quanta in the field 
and the two-level system is in its excited state. In the n- 
excitation subspace, the eigenstates and the eigenvalues 
of the JC Hamiltonian are well known [Hil !T~il | 


6 6 

I ^ +) ) = COS - I n, 0) - sin - | n - 1,1), 

6 0 

I u { ~ ] ) = cos - I n - 1,1)+ sin - | n, 0), 

E i ±} = (n-^jhw± A/2, (3) 

where 



A = (S 2 +4 | g s | 2 n) 1/2 , (4) 


with 6 = h(ui o — ui) the detuning in energy between 
the two-level system and the field mode. The dressed 
states | u^) are entangled atom-field states, and the 
interaction removes the degeneracy present at resonance 
(w = id o) in each subspace with n/0 (the ground state 
is not degenerate and it is not shifted in energy by the 
interaction). Other analogous models frequently used to 
deal with a discrete system interacting with field modes, 
and that can be exactly diagonalized by Bogoliubov-type 
transformations, include an harmonic oscillator interact¬ 
ing with a set of independent harmonic oscillators (see 
for example [2l|). 

The relation between bare and dressed states of the 
JC model can be also expressed in terms of the unitary 
operator [14, 22] 


T = exp 9 (4JV) 1 ^ 2 (a&^ — a^6)| , (5) 

where 8 is an operator defined as 

sind = — (4 | g s | 2 N) 1 ^ A -1 , cosd = — jA _1 , (6) 

with 

A = (<5 2 + 4 | | 2 N) 1/2 . (7) 

Of course, if <5 ^ 0, the operator S 2 + 41< 7 S | 2 AC is self- 
adjoint and strictly positive. Then it admits a positive 
square root, A, which is surely invertible. On the other 
hand, strictly speaking, A/” -1 / 2 is not well defined, since 
zero is one of the eigenvalues of N. However, the inverse 
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of A f (and of its square root) does exist in any subspace of 
T-L which does not contain the vacuum of a and b , which 
are the only subspaces relevant in our analysis. 

The Hamiltonian H , if expressed in terms of the trans¬ 
formed operators, assumes a diagonal form, 

H = huja^a + (fouj — (b^b — , (8) 

where O = TOT~ x is the transformed operator of O (A f 
is invariant under the transformation T: J\f = TAAT -1 = 
A f). Then the dressed eigenstates © can be obtained 
from the corresponding bare ones by application of the 
unitary operator T. Later on in this paper we shall ex¬ 
tend this transformation to the case of a non-Hermitian 
Jaynes-Cummings Hamiltonian model, specifically when 
the coupling constant g s in m is assumed complex. 
This is a situation occurring in the renormalized Lee- 
Friedrichs model in the so-called ghost regime, that is 
when the bare coupling constant is larger than a critical 
value d, [q!, [23| |24|. The Lee-Friedrichs model (23| has 
some similarity with the Jaynes-Cummings model we are 
considering, main difference being that the Lee-Friedrichs 
model contains an infinite (continuous) set of field modes, 
and it is usually restricted to the one-excitation subspace 
only. In the case of the Lee model with imaginary cou¬ 
pling constant (ghost regime), it has been shown that it 
still has a VT (parity-time) symmetry d, Q, and that 
it can be usefully used to describe the light transport 
in an infinite waveguide lattice with a time-dependent 
nonlinearity mu. Also, it has been shown that the 
Hamiltonian of quantum electrodynamics becomes non- 
Hermitian when the unrenormalized electric charge is 
taken to be imaginary, but the Hamiltonian is VT sym¬ 
metric, and the time evolution unitary, if the quadripo- 
tential transforms as a pseudovector [261 ]. 

Following similar lines of the discussion in [8| for the 
Lee model, we can show that also the JC Hamiltonian Jl|) 
is VT symmetric even if the coupling constant g s = ig 
becomes imaginary (and thus the Hamiltonian is not 
Hermitian), provided the real quantity g is such that 
VTgVT = —g . In the next section, we will show that 
the non-Hermitian JC Hamiltonian proves very useful to 
describe the properties of a two-level system interacting 
with a field cavity mode, when a periodic modulation 
of the frequency of the two-level system or of the cav¬ 
ity mode is given to the system. In the first case such 
a modulation can be experimentally obtained by Stark 
shift through the interaction of a two-level atom with an 
external periodic electric field; in the second case it can 
be obtained through a mechanical oscillation of a cavity 
wall, which could be also obtained by a dynamical mir¬ 
ror, that is a wall whose dielectric properties are rapidly 
changed with time, as in experiments recently proposed 
to detect the dynamical Casimir effect [33, dll • 

We now consider the case of an open JC system, in 
which either the transition frequency of the two-level sys¬ 
tem or the frequency of the cavity field mode periodically 
change with time. As mentioned above, the first case 


could be obtained by subjecting the atom to an external 
laser held which periodically modulates its energy lev¬ 
els by time-dependent Stark shift, while the second case 
could be obtained by periodically moving one of the walls 
of a one-dimensional cavity. Both cases allow energy gain 
or loss for the coupled system. In the first case, the JC 
Hamiltonian becomes 

H = b — —^ + Huja l a + g ( ab ^ + <A&) , (9) 

while in the second one we have 

H = Hu ;o b — —^ + hjj(t)a^a + g (atf + a^b) , (10) 

where ui(t) and u>o(t) are prescribed functions of time 
that will be specified later on. We assume that the cou¬ 
pling constant g is real, and thus the time-dependent 
Hamiltonians m and m are Hermitian operators. A 
non-adiabatic modulation of interacting atom-held pa¬ 
rameters has been also recently recognized as a way 
to obtain new phenomena in quantum electrodynam¬ 
ics (QED) involving photon-polaritons coupling |29[, dy¬ 
namical Casimir-Polder forces jjil . [3li ] and generation of 
quantum vacuum radiation [32l . l33l ]. Cases somehow anal¬ 
ogous to the present one have been recently considered 
for the Lee model @, thus restricted to the one-excitation 
subspace only. In our case, they are generalized to an ar¬ 
bitrary n-excitation subspace, including the strong cou¬ 
pling regime near resonance, although for a single cavity 
mode. 


III. SOLUTION OF THE MODULATED 
JAYNES-CUMMINGS MODEL 


We now write the equations of motion for the state am¬ 
plitudes in a generic n-excitation subspace for our mod¬ 
ulated Jaynes-Cummings model, for both cases of peri¬ 
odically driven atomic transition frequency or held-mode 
frequency, described by Hamiltonians m and © re¬ 
spectively. In both cases, in the n-excitations subspace, 
we can write the state in the following form 

I = c'(t) | n - 1,1) + d!{t) | n, 0), (11) 


where, as we said before, the hrst element in the states 
refers to the bosonic held and the second one to the 
fermionic held. Substituting (HU) into the Schrodinger 
equation, we obtain a set of differential equations for the 
coefficients c'(t ) and d'(t). 

In the case of a time-varying atomic transition fre¬ 
quency, using the Hamiltonian ©, we obtain 


- 


‘t/v = 


U} 0 (t) 


+ u(n — 1) 


~c’(t) + 9 -^d\i) (12) 


u 0 (t) 


d'(t) + ^c\t). (13) 
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We now assume that the frequency uio(t) oscillates pe¬ 
riodically in time according to 

uj 0 (t) = uj 0 [1 + 0 cos(Ht)] , (14) 

where ujq is the average frequency, /3 is a complex pa¬ 
rameter and Q is the frequency of the external modula¬ 
tion. From a physical point of view, a complex /3 can be 
thought as a gain or loss term for the atomic energy, as 
we will discuss in more detail later on. 

We now define 

c'(f) = c{t), (15) 

d'(t) = e-^-^ d{t). (16) 

This is equivalent to using the interaction representation 
with an atomic Hamiltonian where the average atomic 
frequency wo appears. Equations (fT2l) and (fT3l) thus yield 

ij t c(t) = ^cos.(fit)c(i) + ^^""^(i), (17) 

i^-d(t) = —^^cos (Q,t)d(t) + e~ i ^°~‘ J - ) *c(f)18) 

at 2 n 

We also define 

c(t) = c(f)e _i ^ sin(nt) . (19) 

Using definitions csd, in and m, equations ED 
and m f° r the amplitudes give 

i—c(t) = 5V5 e ^o-«)t e i^sm(nt) rf ( f ) ( 20 ) 

at a 

i^d[t) = -^cos(nt)d(t) 

_|_ £^ e -*(wo-«)i e -i^sin(nt) c ^\ ^21) 

h 

In the absence of external modulation of the atomic 
frequency, that is j3 = 0, and near resonance uio — the 
solution of the equations above is an oscillation of the 
amplitudes at the Rabi frequency FIr = g^/n/H. When 
the modulation given by (1141) is acting on the system with 
f l Mr, a reasonable assumption because the Rabi fre¬ 
quency is usually much smaller than typical atomic fre¬ 
quencies dH, the amplitudes c(t) and d(t) are expected 
to evolve much more slowly than changes given by the 
modulation frequency f l. In this case we can average 
Equations ( 1 2 () I) and ED over a period relative to the 
modulation frequency H, similarly as it is done in Q, 
bringing slowly-varying terms outside the time integrals. 
Using 


the averaged equations at resonance wq — w are thus 



As a specific case, if we put n — 1 in our results we 
recover previous results in [9( for the one-excitation sub¬ 
space (Lee-Friedrichs model). 

We can now compare Eqs. (1231) and (l24l) with the anal¬ 
ogous equations obtained in the static case (/? = 0), that 
is a Jaynes-Cummings model at resonance with time- 
independent atomic and field frequencies, described by 
the Hamiltonian dT|) 

4 c(f) = <25) 

4 d(f) = <26) 

If we choose the complex parameter /3 in such a way 
that 

■ 4 -^)=., ( 27 ) 

equations (l23l) and (l24l) for the modulated case coincide 
with those for the static case given by Eqs. (l25l) and (l26l) 
with 

g s = ig € C , (28) 

that is with an imaginary coupling constant making non- 
Hermitian the (time-independent) Hamiltonian (HD. We 
point out that equation (1271) has a complex solution given 
by — —2.14 + 1.42 i, yielding a complex value of /3. 

Same considerations apply to the case of a modulated 
field frequency u(t), described by the Hamiltonian (fTTTT) . 
Also in this case we may assume a periodic modulation 
of the frequency of the cavity mode, 

uj(t) = U) (1 +/3cosfli) , (29) 

with /3 a complex parameter. Using ED , we can write 
the equations for the amplitudes and use definitions (fTTH) 
and ED. Defining now 

c(t) = c(t)e~ iLsin ^ , (30) 

we obtain a set of differential equations for the ampli¬ 
tudes analogous to (l23l) and (l24l) . Thus, also in this case, 
comparison with Eqs. (l25l) and (l26l) shows that by the an 
appropriate choice of the system parameters such that 


1 

2-7T 


0 —ir sin x 


dx = J 0 (T). 


( 22 ) 



(31) 


where T is a complex number and Jq(z ) is a Bessel func- we can simulate, also in this case, the Hermitian JC 
tion [35j], and observing that Jq(z) is an even function, Hamiltonian with a modulated cavity frequency by a 
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static non-Hermitian JC Hamiltonian with an imaginary 
coupling constant. 

A complex value of the parameter [3 in m as obtained 
by the solutions of equations (12711 and (IdTl) . could result 
from several microscopic physical mechanisms involving 
energy exchange from or to the system, for example from 
the interaction of the atom with an external environment 
whose properties are modulated in time. Such a mod¬ 
ulated environment could be a set of infinite harmonic 
oscillators with a continuous energy spectrum in out of 
equilibrium conditions (that can exchange energy with 
the atom), by a modulated (i.e. oscillating) optical cav¬ 
ity or by a dynamical photonics crystal. 

These results thus give a physical meaning to the non- 
Hermitian Jaynes-Cunnnings Hamiltonian © with an 
imaginary coupling constant g s = ig , g being a real 
quantity. In other words, our results clearly show that 
a Jaynes-Cummings model with an imaginary coupling 
constant, yielding a non-Hermitian Hamiltonian, can be 
a useful model to simulate a realistic Hermitian Jaynes- 
Cummings system with some external modulation of the 
physical parameters, specifically the cavity mode fre¬ 
quency or the atomic transition frequency, at least when 
gain-loss effects want to be efficiently described. This 
gives a physical support to further investigate a non- 
Hermitian Jaynes-Cummings model, from both physical 
and mathematical aspects. 

In recent years, it has been discussed that physical sys¬ 
tems described by non-Hermitian VT symmetric Hamil¬ 
tonians can be investigated in the framework of the more 
general formalism of pseudo-Hermitian operators [l8l . l36j . 
As we shall discuss in the next Section, pseudo-operators 
are defined by introducing a particular deformed version 
of the canonical commutation relations. The meaning of 
the pseudo-operators and their relation with VT sym¬ 
metry has been recently considered, both from a purely 
mathematical point of view and also in view of possible 
applications to the so-called pseudo-Hermitian quantum 
mechanics, [(37). Pseudo-Hermitian quantum mechanics 
has received a lot of interest since the discovery of the 
fact that the well known cubic Hamiltonian H = p 2 +ix 3 , 
has purely real and discrete spectrum ,6). The interest 
on this subject is also motivated by the fact that the 
spectral properties of VT symmetric Hamiltonians fol¬ 
low from their pseudo-Hermiticity [H;|. 

The relevance of pseudo-Hermitian quantum mechan¬ 
ics in the framework of quantum optics and open quan¬ 
tum systems has been discussed. For example, it has 
been recently shown that non-Hermitian interaction be¬ 
tween atoms and the electromagnetic field can be ap¬ 
propriately investigated by exploiting the methods of 
pseudo-Hermitian quantum mechanics [36| . 

In the next section, on the basis of the previous re¬ 
sults showing the physical interest of the non-Hermitian 
Jaynes-Cummings Hamiltonian with a purely imaginary 
coupling constant, we will discuss a deformed version of 
the Jaynes-Cummings Hamiltonian, in which the bosonic 
and fermionic modes are replaced by their pseudo- 


bosonic and pseudo-fermionic versions, [38[ . This inter¬ 
est is motivated by several recent results, see for instance 
fTol . [39l l4pj | and [HI for a recent review, where many 
systems introduced along the years in connection with 
VT or pseudo-symmetric Hamiltonians, where complex¬ 
valued parameters are introduced, have been shown to 
be expressible in terms of pseudo-bosonic or pseudo- 
fermionic operators. In doing so, the extra bonus one 
gets is that the eigensystems of these Hamiltonians, and 
of their adjoint, can be easily deduced, following a some¬ 
how standard procedure, procedure which will be used 
to diagonalize the Hamiltonian in (15^1) below, generaliz¬ 
ing to the case of pseudo-operators the procedure used in 
the Hermitian case. This will have also very interesting 
mathematical consequences, as we will see, since pseudo- 
bosonic and pseudo-fermionic structures are mixed along 
the way. Then we conclude that the deformed model 
of the Jaynes-Cummings Hamiltonian we introduce be¬ 
low is thus physically well-motivated and permits one to 
appropriately investigate the non-Hermitian interaction 
between a two-level system and the electromagnetic field 
in a cavity. 


IV. THE NON-HERMITIAN 
JAYNES-CUMMINGS HAMILTONIAN 

Pseudo-bosons and pseudo-fermions are defined per¬ 
forming a particular deformation of the canonical com¬ 
mutation relations for bosons and fermions. In partic¬ 
ular, pseudo-bosons arise from the canonical commuta¬ 
tion relation [a, ad] = 1 after replacing a' by another 
operator A satisfying the relation [a, A] = 1, with A in 
general different from a'. More recently, a similar proce¬ 
dure to define pseudo-fermions has been carried out, and 
their properties have been analyzed again both from a 
mathematical and from a physical point of view [H] • We 
now introduce our model for a non-Hermitian Jaynes- 
Cummings Hamiltonian expressed in terms of pseudo- 
bosonic and pseudo-fermionic operators, and we discuss 
its diagonalization. 


A. The deformed Jaynes-Cummings Hamiltonian 

The Hamiltonian of our model, which can be consid¬ 
ered as the extension of the Hamiltonian CD considered 
in Section [Til can be written as 

H a = tujj 0 ^ 'C a c a — + fojjD a d a + ed a C a + e*D a c a , 

(32) 

acting on a Hilbert space 7i := 'Hb®'Hf, where Hf = 
C 2 (fermionic sector) while TLb is infinite dimensional 
(bosonic sector), e indicates the coupling constant (we 
are using different symbols in order to avoid confusion 
between the Hermitian and non-Hermitian cases). The 
suffix a is a parameter which, when is sent to zero (in 
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an appropriate way) produce a self-adjoint Hamiltonian 
H 0 = Hq. On the other hand, since for a ^ 0 we have, in 
general ct ^ C a and <fi a ^ D a , we expect that H a ^ H 1. 
More in details, following a general result discussed in 
[38j . what we have in mind is that c a , C a , d a and D a 
are (at least formally) similar to two pairs of bosonic (c 
and ct) and fermionic ( d and d') operators, and that this 
similarity is implemented by some invertible operator Q a , 
possibly unbounded with unbounded inverse, such that 
Q a converges to the identity operator (in some suitable 
topology) when a goes to zero. We will briefly return to 
this point in Section liVBI 

We assume here the following rules 

[d a 0 1/,-Da 0 1/] = Ifc 0 1/ =: 1 , 

{l b 0 c a , l b 0 C a } = 11, (33) 

while all the other commutators are zero. Hence the pair 
( d ai D a ) satisfies the pseudo-bosonic commutation rules, 
while ( c ai C a ) behaves as pseudo-fermions. 

With a simple extension of the procedure discussed in 
Section |TT] (see also 0 ), we can rewrite H a in a diagonal 
form. For that we first introduce a global non self-adjoint 
number operator, analogous to the excitation-number op¬ 
erator ([2]), 

N a = D a d a + C a c a , (34) 

and a map T a defined as follows: 


[38i |. We first assume that two non-zero vectors Ao and 
ipo do exist in T~Lb such that, if r)o an d Ao are two vectors 
of fermionic Hilbert space T~i / annihilated respectively by 
c a and C£, we have 


da 0 1/^ $0,0 = 

(lb 0 C a ) $o,0 = 0, 

(37) 

as well as 



(d+ 0 1/) $o,o = 

(lb 0 j $ 0,0 = o, 

(38) 


where $ 0 ,o : = An 0 Vo and $ 0 ,o := 'Ao 0 Ao- 

Following |38j, it is convenient to assume that Ao and 
Ao belong to a dense domain V of Tib, which is stable 
under the action of d a , D a , and their adjoint. Of course, 
nothing has to be required to the fermionic operators and 
to the related vectors, since T~Lf is a finite dimensional 
vector space. 

If such a T> exists, then we can use the two vacua $o,o 
and d'o.o to construct two different set of vectors, := 
{$n,fc, n > 0, k = 0,1} and F^, := n > 0, k = 

0,1}, all belonging to V <g> 'Hf , as follows: 



T a = exp |-6» a (4|e| 2 A r a ) 1/2 {ed a C a - e*-D a c a ) j , 

(35) 

where 9 a is an operator defined, in analogy with ([6]), via 
the conditions sind a = —(4|e| 2 A r a ) 1 / 2 A“ 1 and cos 0 a = 
—(5A" 1 , where S = h{ujQ — uj) is the detuning between the 

energies of the two fields, and A a = (<5 2 + 4|e| 2 -/V ct ) , 

which is clearly invertible, apart from when restricted to 
the zero-excitation subspace (we are not however inter¬ 
ested in this subspace, because our Hamiltonian is al¬ 
ready diagonal in this subspace). 

It is easy to see that the new operators d a 0 1/ = 
T a (d a ®TLf)T- 1 , D a ® 1/ = T a {D a ® l/)! 1 ” 1 , U0c a = 
T a (l b 0 CajT- 1 and i b 0 C a = T a { l b 0 CA)! 1 " 1 still 
obey the same rules as the operators without the hat: 
hence they are tensor products of pseudo-bosonic and 
pseudo-fermionic operators. Most important, in terms of 
these operators the Hamiltonian H a turns out to be in a 
diagonal form 

H a = (twj - A a ) ^ C a c a - 0 + huDada , (36) 

/ . \l/2 

where A a = ( S 2 + 4\e\ 2 N a J = A a , since N a = N a . 

The eigenvectors of H a (and of H' a ) can now easily 
computed adapting to the present situation the general 
framework of deformed canonical commutation relations 
and canonical anti-commutation relations discussed in 


and 


4' n,k = 


Vn\ 

1 

— i=t 

\fn\ 


cl ) $ 


- 0,0 


dl Ao ) 0 ( ct Q Ao =: An 0 Afe ,(40) 


with obvious notation, where n = 0,1,2,... and k = 0,1. 
It is now easy to check that 

H a *&n,k = En,A;$n,/c; ■ (41) 

where 


E„ i- — 


Huj — (d 2 + 4|e| 2 (n + k + 1)) 


1/2 ^k-^J+hun. 

(42) 

Incidentally, we see that E n ^ is surely bounded from 
below in n, independently of the choice of the parameters. 
As for k = 1, we see that this is true only if Hlo > 2|e|, 
which is therefore to be assumed to give our Hamiltonian 
a physical meaning. 

Equation (Hill ensures us that F^ and F^, are biorthog- 
onal. Indeed, if the normalizations of $o,o and $o,o are 
chosen in such a way that ( < I > o,o, $o,o) = 1, then 


$n,fc, $m,Z j — \ Tni {Vkj A• 

(43) 
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Here (., .) Ub and (., .) Hf are respectively the scalar prod¬ 
ucts in 'Hb and in 'Hf. 

A serious problem, as it happens quite often in this 
kind of problems, is whether E^ and E^, are bases of T~i 
or not. Or, at least, if they are (J-quasi bases, see below, 
for some Q dense in 'H. In general, the answer is hard to 
find. However, there is at least one situation where this 
can be done, as we are going to discuss next. 

B. What if H a is similar to a self-adjoint operator 

We consider now the following self-adjoint Hamiltonian 

H 0 = huj 0 + Hujd^d + edc 1 + e*d)c , (44) 

where [d, d 1 '] = lb, {c, c’ 1 '} = 1/, c 2 = 0, and all the other 
commutators are zero. It is clear that Hq = Hq. It is 
also evident that this can be obtained from H a in (l32l) re¬ 
placing pseudo-bosonic (fermionic) operators with their 
standard counterparts or, in view of our previous com¬ 
ment, sending a to zero. Hq is the Jaynes-Cummings 
Hamiltonian essentially discussed in GiQ , and it can 
be diagonalized as we did for H a . In particular, intro¬ 
ducing new bosonic and fermionic operators by means 
of the unitary operator To, see (l35l) . the eigenstates of 
Hq can be easily deduced. More explicitly, introducing 
d ® 1/ = T 0 (d ® 1/)T 0 -1 , and R 0 c = T 0 (lb ® c)T 0 -1 , 
the Hamiltonian can be rewritten as 

Hq = (juu — ^ + huicftd , (45) 

~ „ - \V 2 

where A 0 = f S 2 + 4|e| 2 -/V 0 J = A 0 , since N a = N a . 

The eigenstates <f>° k of Hq can be defined easily: first 
we introduce the vacua of d and c, ipQ £ Hb and t)q £ Hf: 
dip o = ct)q = 0. These vectors surely exist, as it is clear. 
Then we construct, more solito, = -^= (E) n ifQ and 
f/° = c^q, and finally we define <f>° k = (p° n <g> r) k , n > 0 

and k = 0,1. Since iJo<£° fc = E nik ®° >k , H 0 , H a and H\ 
are all isospectral. Hence intertwining operators between 
them are expected to exist, |4lj| . For this reason, it is not 
a strong assumption to assume here that a self-adjoint 
operator S a = SJ t does exist, at least on some dense 
subset of H, H a S a = S a HQ. Then, if we further assume 
that S a is invertible, we see that H a = SuHqS^ 1 . Of 
course we also have H^ a = S~ 1 HqS ol and, more important 
for us, 

$n,k = s a ^ k , 9 n>k = S,, ( 46 ) 

at least if e D{S a ) fl D(S~ 1 ) =: Q, where D(S a ) is 
a dense domain where the operator S a is defined, and if 
the multiplicity of each E n)k is one. These equations are 
in agreement with formula (H51) : this is because the set 


= {$° k } is an orthonormal basis for H. But we get 
more than this: in fact, if S a and S are both bounded, 
then Tjj, and J-^ are biorthogonal Riesz bases, which is 
the best we can have (from a technical point of view) 
when we lose orthonormality, [42l-l44|. However, when 
S a or S^ 1 , or both, are not bounded, the situation is 
not so simple and, in fact, only some weak resolutions of 
the identity can be deduced [38|. In particular, assuming 
that Q is dense in H, we can conclude that Tk and 
are biorthogonal Cl-quasi bases. We refer to [38|| for more 
information about the properties of these sets. Indeed, 
they give rise to several interesting properties which have 
been investigated in recent years. Here we just want to 
mention that, in the present context, when we say that 
E j, and E^ are C/-quasi bases we mean that, taken two 
arbitrary vectors / and g in Q, the following equalities 
are satisfied: 

if,g) = J2 (/> ^,k) (^n,k,g) = J2 (/. 4V*} (*",*>$) • 

n,k n,k 

(47) 

From (l46l) we also deduce that 

= Sl^ nt k, (48) 

so that = S~ 2 ^n,k- These equations show that 

(i) a metric operator, see [38| can be introduced in the 
game, at least formally, which is positive, self-adjoint, 
invertible, and maps one set of eigenstates into the other; 

(ii) the following equalities hold, on some domain 

= £!*».* ><*».*!’ 5- 2 = ^|^„, fe )('hn, fc |. (49) 

n,k n,k 

Here (|/ >< g\) h :=< g,h> /, for all f,g,h £ H. More 
mathematical details on the general framework are dis¬ 
cussed in jHI], to which we refer the interested reader. 

From a physical point of view, the advantage of the 
diagonalization of the non-Hermitian Jaynes-Cummings 
Hamiltonian obtained in this section through the opera¬ 
tor (l35l) . in terms of pseudo-bosonic and pseudo-fermionic 
operators, is that it allows us to have an explicit form of 
the diagonal Hamiltonian and of the dressed states, that 
can be obtained from the bare ones by applying the oper¬ 
ator (|J5|) . in analogy with the usual Hermitian case. This 
can be helpful in investigating many physical systems 
that can be simulated by a non-Hermitian Hamiltonian, 
such as the system considered in Section HTT] , modulated 
optical lattices and the Lee model in the ghost regime 

V. CONCLUSIONS 

In this paper we have considered two time-dependent 
Jaynes-Cummings-type Hamiltonian models, describing 
a two-level system interacting with a single cavity mode. 
More in details, we have considered a periodically driven 


transition frequency of the two-level system and a peri¬ 
odic change of the frequency of the cavity mode. These 
two interesting situations could be obtained, respectively, 
by a dynamic Stark shift of the energy levels of an atom 
by a laser field, and by a motion of the conducting wall 
of a cavity or even a dynamical mirror. We have first 
shown that, for appropriate choices of the system pa¬ 
rameters, both cases can be simulated with a static non- 
Hermitian VT symmetric JC Hamiltonian with an imag¬ 
inary coupling constant. This result suggests the rele¬ 
vance, even from a physical point of view, to consider the 
non-Hermitian JC model. With this in mind, we have an¬ 
alyzed our system by exploiting a more general procedure 
for treating non-Hermitian Hamiltonians. We have in¬ 
troduced a model of deformed Jaynes-Cummings Hamil¬ 
tonian expressed in term of pseudo-operators, and ex¬ 
tended the well-known unitary operator diagonalization 
of the Hermitian JC Hamiltonian with a invertible opera¬ 
tor to the non-Hermitian case in terms of pseudo-bosonic 
and pseudo-fermionic operators. This has allowed us to 


obtain explicitly its diagonal form and its eigenstates 
and eigenvalues, analogously to the known Hermitian 
case. This has also very important consequences from 
the mathematical point of view, since pseudo-bosonic 
and pseudo-fermionic structures are shuffled along the 
way. Mathematical and physical implications and appli¬ 
cations of our results have been also discussed in detail. 
The deformed JC model introduced could be used to fur¬ 
ther investigate the interaction between atomic systems 
and the electromagnetic field, including damping or am¬ 
plifying processes, which is of fundamental importance 
for example in quantum optics. 
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